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ABSTRACT 

An algorithm to forecast occurrence of four life-stage categories of the alfalfa weevil, Hypera 
postica (Gyllenhal), was derived from data collected in fields of seed alfalfa, Medicago sativa 
(L.) in southern Alberta. The algorithm assumes a linear developmental response to mean daily 
temperatures above a threshold of 10°C. Overwintering adults were active after the accumula¬ 
tion of 100 degree-days above 10°C (DDio), and were scarce by 250-300 DDio- Early larvae 
(instars 1 + 2) were found beginning at 120 DDjo and their numbers peaked at 200 DDio. Late 
larvae (instars 3 + 4) were present beginning at 160 DDjo and their numbers peaked at 350 
DDio New generation adults appeared after 500 DDio. In southern Alberta, alfalfa seed pro¬ 
duction is frequently combined with honey production. This algorithm enables producers to 
forecast the occurrence of the most damaging stage of alfalfa weevils which may require con¬ 
trol with insecticides; the advance notice enables optimal timing of treatment and also allows 
apiarists to minimize pesticide mortality by moving or confining their bees. 

Key words: Thermal units, degree days, simulation, phenology 

INTRODUCTION 

Pest control measures are economically justifiable only if their benefits exceed their cost 
(Stem et al. 1959). Normally, cost is the sum of pesticide purchase plus its application, and 
benefit is measured by reduced yield loss. However, other factors may enter the cost:benefit 
equation; this occurs in alfalfa production in southern Alberta, where seed producers frequently 
obtain additional income by charging apiarists to place honeybees (Apis mellifera L. (Hy- 
menoptera: Apidae)) in their fields. Thus, strategies to control pest insects in alfalfa seed fields 
must acknowledge the susceptibility of honeybees to many pesticides. 

One approach to reconciling alfalfa pest management with apiculture is to give apiarists suf¬ 
ficient warning to cover the hives, or to move them, before pesticides are applied. This approach 
requires a method of forecasting the occurrence of the pest population. This paper presents a 
simple method of doing so, which is based on observed correlations between phenological 
events and degree-day accumulations. 

In southern Alberta, the alfalfa weevil is a serious insect pest of alfalfa, feeding on shoots, 
flower buds and foliage during prebloom to early bloom (Hamlin et al 1949). If not controlled 
it can severely reduce seed yield. Adult alfalfa weevils spend the winter in protected locations 
in alfalfa fields or in litter nearby. Overwintered adults become active about the time the first alf¬ 
alfa shoots appear in the spring. They feed for a few days, mate, and begin ovipositing. Peak egg 
density usually occurs in late May or early June, but eggs can be found during most of the sum¬ 
mer. 

Rates of pre-imaginal development are temperature-dependent. Egg incubation takes 4 to 21 
days. Larvae develop through the four instars in about 3 to 4 weeks. Instars 1 and 2 (early lar¬ 
vae) feed within the tightly-curled developing leaves and buds; instars 3 and 4 (late larvae) feed 
on expanded leaves. 

Feeding damage is most obvious in mid- to late-June, concurrent with peak densities of late 
larvae. Defoliation is most severe toward the terminals. This type of feeding results in loss of 
foliage, flower buds, and nutrients, with consequent delays in plant growth and development. 
Quantity and quality of the seed yield may be reduced. 

Late larvae inflict the greatest damage, so agricultural losses can be expected to start accu¬ 
mulating at the onset of this life stage (Dennis et al 1986). Consequently, the timing of control 
tactics for the alfalfa weevil should be optimized, based on an understanding of physiological 
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Figure 1 . Observed and predicted proportion of alfalfa weevil in relation to DDjo. a) Over¬ 
wintering adults (life stage = 1): o and solid line. New generation adults (life stage = 4) A and 
dotted line; b) first and second instars (early larvae; life stage = 2); c) third and fourth instars (late 
larvae; life stage = 3). Lines are predictions made using Equations 1-3 in text. 
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and behavioral processes (Harcourt 1981; Whitford and Quisenberry 1990), to target this life 
stage. This optimization requires information on the seasonal abundance and time of occurrence 
of the immature stages of the alfalfa weevil through a combination of population monitoring 
and phenological models (Schaber and Richards 1979). 

In southern Alberta, chemical insecticides are the primary method of alfalfa weevil control. 
Efficient use of insecticides requires that applications should be timed to target the first late lar¬ 
vae, i.e. after they have moved to exposed positions on the leaves, but before they have caused 
much damage. This timing requirement establishes a potential conflict between pest control 
and honey production. A reliable method of forecasting the occurrence of late larvae could re¬ 
solve this conflict. 

This study was conducted to develop a technique that would predict the appearance of late 
instar alfalfa weevil larvae in seed alfalfa fields, using southern Alberta field data. This tech¬ 
nique would enable better timing of insecticide applications in relation to insect development, 
and allow seed producers and apiarists more lead time to protect pollinators. 

The temperature-dependence of insect development dictates that developmental models be 
based on thermal-unit accumulation. Simple methods for modelling phenology based on field 
data are available, and can be used to develop realistic models in the absence of detailed data on 
insect development processes (Hudes and Shoemaker 1988; Kemp and Onsager 1986; Kemp 
et al 1986; Ly syk 1989). Degree-day accumulation has been used to predict peak hatch and sub¬ 
sequent activity of alfalfa weevil in forage alfalfa in southern Ontario (Harcourt 1981). How¬ 
ever, Tauber et al (1988) have suggested that the phenology of an insect species can vary among 
geographic regions due to adaptation of thermal biology to local climatic conditions, so another 
objective was to compare phenology of the alfalfa weevil populations in southern Ontario and 
southern Alberta. 


METHODS 


Algorithm development 

The algorithm was developed using phenology data obtained from four research plots at the 
Agriculture and Agri-Food Canada Research Centre (AACRC) at Lethbridge, Alberta. Alfalfa 
weevil abundance was determined by taking five sweeps per plot with a 38-cm net (Johansen 
et al. 1979). Each plot was sampled every one to three days from 3 June to 12 August 1985; 
12 May to 8 August 1987; 3 June to 20 July 1988; 3 June to 4 July 1989, and was sampled 
weekly from 4 June to 23 July 1990. Daily maximum and minimum temperatures (°C) were ob¬ 
tained from the AACRC meteorological station. 

Alfalfa weevil phenology was modeled by correlating phenological events with degree-day 
accumulations using a developmental threshold of 10°C. This threshold was used because, al¬ 
though alfalfa weevil eggs hatch at 8°C, the resulting larvae do not survive even if subsequently 
exposed to a higher temperature (Guppy and Mukerji, 1974). Accumulated degree-days above 
10°C (DDio) from 1 January were calculated for each year by sine-wave integration (Allen 
1976). Accumulated DDio on each date was rounded to the nearest 10, and samples from the 4 
plots were grouped according to these rounded values. The proportion of alfalfa weevils which 
were adults, early larvae and late larvae were calculated for all such grouped samples, and these 
proportions were related to the rounded DDio accumulations by non-linear regression (Proc 
NLIN, SAS Institute 1989) as outlined below. 

To provide a standardized estimate of when the weevils and larvae were most abundant, and 
to establish correlations between phenological events and DDio accumulations, their relative 
abundance was calculated for each plot on each day in each growing season by summing the 
number of weevils collected and dividing by the greatest number collected in one day for that 
plot in that year. These relative abundances were grouped across plots and years by the rounded 
degree-day values. The correlations between relative abundance and DDjo accumulations were 
used to develop an algorithm to predict the appearance of the different life stages. 

Alfalfa weevil phenology was divided into four stages (i): overwintered adults (i = 1), i.e., 


Proportion of population in life stage i or greater 
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Figure 2. Proportion of alfalfa weevil in a) life stage 2 or greater, b) life stage 3 or greater, and 
c) life stage 4. Solid lines are model predictions using equation 2 and parameters estimates 
given in Table 1. 
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those present before the larval peak at 350 DDio; early larvae (i = 2); late larvae (i = 3); and new 
generation adults (i = 4), i.e, those occurring after the larval peak at 350 DDio. 

The algorithm was developed as outlined below (Hudes and Shoemaker 1988). For each 
rounded DDio value, the following calculations were made: 

F 2 = (n 2 + n 3 + / 2n 

F 3 = (n 3 + / 2n 

F 4 = n 4 /2n 

where ni(i=2-4) is the number of weevils in life stage i, 2n is the total number of weevils, and 
Fj is the proportion of insects in life stage i or later. Note that because all insects are in a life stage 
equal or greater than the overwintered adult stage, Fi = 1. 

The time trends in each Fj were modelled using equation 1. 

F, = [l (1) 

Non-linear regression was used to obtain estimates of the parameters (aj, bj; i = 2 - 4). The 
variable tj is a scaled estimate of thermal time calculated for each life stage as: 

_ DDio - DD m i n (,-) 

ij —- ( 2 ) 

In equation 2, DD min o) and DD max (|) represent the approximate value of DDio for the be¬ 
ginning and end of life stage i, obtained by inspection of die data. 

The functions describing time-change in \ were then used to predict the proportion of insects 
in each life stage (pi): 

Pi = l-F 2 , 

P2 = F 2 - f 3> ( 3 ) 

P3 = F 3 -F 4 , 

P4 = F 4 . 

Algorithm validation 

Independently-obtained data were used to validate the algorithm. These were collected from 
plots in Brooks, Rosemary and Rolling Hills, Alberta (ca. 130 km N.E. of Lethbridge) by 
WestAg, Inc., a pest management scouting company. These plots were sampled weekly for 
up to 13 weeks in 1984-1988, starting the last week of May and continuing through August. 
Samples in a specific plot in each year were taken at approximately the same time of day to min¬ 
imize any effects of insect diurnal cycle on sampling efficiency (Johansen etal. 1979). The data 
consisted of weekly mean numbers of adults and of early and late larvae. The numbers of fields 
sampled were 48, 40, 42, 48 and 42 for 1984 through 1988. Sample counts were tabulated 
weekly. Daily maximum and minimum temperatures were obtained for each year from the 
Alberta Special Crops and Horticultural Research Center in Brooks, Alberta, and the DDio 
accumulation from January 1 were calculated by sine-wave integration (Allen 1976). The DDio 
accumulation in the middle of each week was matched with the weekly insect data. The aver¬ 
age proportion of weevils in the adult, early larval, and late larval stages was calculated for each 
week and compared graphically to the algorithm predictions. 
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Figure 3. Proportion of alfalfa weevil in commercial seed alfalfa fields. Symbols are observed 
mean proportion ± 1 standard error, and lines are predictions made using Equations 1-3. a) 
adults, solid line is overwintered generation (life stage = 1) and dotted line is new generation 
(life stage = 4); b) first and second instars (early larvae; life stage = 2); c) third and fourth instars 
(late larvae; life stage = 3). 
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RESULTS 


Seasonal phenology of alfalfa weevil 

Adult weevils were first found at about 100 DDm (Fig. la) and at that time were the only 
stage collected. The proportion of adult weevils declined from about 250 DDjo and they had be¬ 
come scarce by 350 DDm. Early larvae appeared at about 120 DDm (Fig. lb), reached a max¬ 
imum at about 200 DDm, and then declined slowly. Late larvae appeared at ca. 160 DDm (Fig. 
lc), peaked near 450 DDm, and then declined. New generation adult weevils began to appear 
after 400 DDm, and increased steadily to nearly 100% of the population by 900 DDm- 

Algorithm output 

Figure la-c illustrates the good agreeement between the algorithm output and the Lethbridge 
data. Estimates of the regression parameters are listed in Table 1. Coefficients of determination 
(r 2 ) were 0.97,0.84, and 0.74 for life stages 2, 3 and 4 respectively. Model predictions of pro¬ 
portions of insects in life stage i or greater are overlaid with observations in Fig. 2a-c 

Algorithm validation 

Observed and predicted values of the proportion of weevils in each life stage for the test data 
are shown in Fig. 3 a-c. The predicted proportions of weevils in the adult life stages (1 and 4) 
seem reasonable (Fig. 3a). For early larvae, the algorithm predicted the timing of the peak 
population (Fig. 3b), and in general captured the seasonal population trends, but overall un¬ 
derestimated the proportion of weevils in this life stage. For late larvae (Fig. 3c) the algorithm 
predicted the start of the stage and captured the essence of the seasonal trends, but tended to 
overestimate the proportion of insects in this stage. 

DISCUSSION 

The phenology we observed differed from that recorded in Ontario (Harcourt 1981). In On¬ 
tario, most of the feeding damage occurred between 260 and 335 Degree-days base 9°C, 
whereas in this study, the damaging stage extended beyond 500 DDm- These differences may 
reflect adaptation of the insect to a different cropping pattern, i.e. seed alfalfa vs forage alfalfa, 
or to a different climate. Harcourt’s (1981) model does not seem to apply to the Alberta popu¬ 
lation of alfalfa weevil. 


Table 1 

Estimated parameters for component equations of the algorithm simulating alfalfa weevil population phenology. 


Equation (1 ) a Equation (2) b 


Life Stage (i) 

a 

b 

r 2 

DDmin 

DDmax 

2 

18.15 

11.01 

0.97 

120 

400 

3 

4.84 

1.00 

0.84 

160 

800 

4 

5.05 

29.81 

0.74 

200 

900 


a Nonlinear regressions of the temporal change in proportions of observed insects in life stage (i), using 

F, = [1 - e 

b Scaled estimates of thermal time for life stage (i), calculated as: 

_ DD\ 0{[) - DD m - m{[) 

ij — - 

DD max (o - DD m j n m 


Details are given in the text. 
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Although qualitatively similar, there is some discrepancy between the predictions and the 
validation data. Several factors can affect the accuracy and reliability of the algorithm. First, 
separating alfalfa weevil larvae into life stages is somewhat subjective, particularly when dif¬ 
ferentiating between second and third instars, and variation among observers can affect the ac¬ 
curacy of sampling for pest insect populations (Shufran and Raney 1989). 

Another possible source of the discrepancies results from the differing ways in which sam¬ 
ples were handled. The alfalfa weevils and larvae captured in the sweep net at Lethbridge were 
placed in paper bags and stored at-40°C until counted, whereas in the validation data they were 
counted in the field immediately after capture. These differences in method, coupled with the 
subjective assignment to stages, could affect the performance of algorithm. 

Another possibility is that the use of degree-days may not be strictly applicable because it dis¬ 
regards both the nonlinearity of the developmental rate function (e.g. Lactin and Holliday 
1993), and the ability of insects to control body temperature behaviourally (Huey and King- 
solver 1989), and can not account for the effects of transient unfavourable conditions (Lactin 
1992). 

Finally, insecticides were used in the commercial fields surveyed by WestAg, Inc., and dif¬ 
ferential mortality among the life stages may acount for the bias between algorithm predictions 
and field data. Early larvae feed mostly within the flower buds and are largely protected from 
contact insecticides, whereas late larvae feed on expanded foliage and are not (Johansen et al. 
1979). These differences in exposure risk could bias the population structure of the validation 
fields, compared to that predicted by the algorithm, which was based on observations from in¬ 
secticide-free fields. 


CONCLUSION 

This paper outlines the development of an algorithm to predict the occurrence of damaging 
stages of the alfalfa weevil. Although there is a consistent bias in the estimate of proportion of 
these stages in the population, the algorithm estimates the timing of stages quite well, and can 
provide sufficient advance warning to optimize the timing of insecticide application, and thus 
allow apiarists to remove or confine their bees. 
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